Discrete Nonlinear Schrodinger Equations Free of the Peierls-Nabarro Potential 
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We derive a class of discrete nonlinear Schrodinger (DNLS) equations for general polynomial 
nonlinearity whose stationary solutions can be found from a reduced two-point algebraic problem. 
It is demonstrated that the derived class of discretizations contains subclasses conserving classical 
norm or a modified norm and classical momentum. These equations are interesting from the physical 
standpoint since they support stationary discrete solitons free of the Peierls-Nabarro potential. As a 
consequence, even in highly-discrete regimes, solitons are not trapped by the lattice and they can be 
accelerated by even weak external fields. Focusing on the cubic nonlinearity we then consider a small 
perturbation around stationary soliton solutions and, solving corresponding eigenvalue problem, we 
(i) demonstrate that solitons are stable; (ii) show that they have two additional zero-frequency 
modes responsible for their effective translational invariance; (iii) derive semi-analytical solutions 
for discrete solitons moving at slow speed. To highlight the unusual properties of solitons in the new 
discrete models we compare them with that of the classical DNLS equation giving several numerical 
examples. 

PACS numbers: 03.40.Kf, 63.20Pw 
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I. INTRODUCTION 

Recent theoretical and experimental results have 
demonstrated that the fundamental properties of waves 
can be precisely engineered by introducing a periodic 
modulation of the medium characteristics. In particular, 
there appear unique possibilities to control the propaga- 
tion of light in photonic structures with a periodic mod- 
ulation of the optical refractive index [1], and manage 
the flow of atomic Bose-Einstein condensates (BEG) in 
periodic potentials [2]. Additional flexibility in tailoring 
wave dynamics can be realized through nonlinear self- 
action, that may appear due to various physical mecha- 
nisms such as light-matter interactions for optical beams 
or atom-atom scattering in BEG. One of the key effects of 
nonlinearity is the suppression of the natural tendencies 
of localized wave packets to broaden due to dispersion 
or diffraction, supporting the formation of optical lattice 
solitons [3, 4] and self-localized atomic states [5, 6]. 

The excitation of lattice solitons and their dynamics 
can be described by the discrete nonlinear Schrodinger 
(DNLS) equation when the wave propagation is primar- 
ily defined by tunneling between neighboring potential 
wells [3, 7]. This special regime of energy flow results in 
special properties of discrete solitons. It was predicted 
theoretically [8-12] and observed experimentally [13, 14] 
that the discrete solitons can freely propagate through 
the lattice below a certain energy threshold, whereas at 
higher energies the solitons become trapped at a partic- 
ular lattice site. This phenomenon is related to the pres- 



ence of the self-induced Peierls-Nabarro potential (PNp) 
barrier, defining the difference of energies between the 
solitons whose centers are positioned at a lattice site or 
in-between the neighboring lattice sites. 

The ability to enhance or suppress PNp may allow for 
precise control over transmission of higher-energy wave 
packets through the lattice. The PNp vanishes for dis- 
crete solitons in the framework of the Ablowitz-Ladik 
(AL) equations [15, 16] which are integrable. One of 
the features of the AL model is nonlocal nonlinearity, as 
the medium response is defined by the wave amplitudes 
at the neighboring lattice sites, in contrast to the local 
on-site nonlinearity that leads to strong PNp [8]. How- 
ever, the AL model has not been directly connected to 
a specific physical application. On the other hand, more 
recently, it was demonstrated that nonlocal nonlinearity 
of a more general type representing real physical systems 
can indeed lead to strong reduction of PNp [17, 18] for 
lattice solitons. It was also shown that PNp can be re- 
duced (and even reversed) in the case of a saturable non- 
linear response [19]. 

In this paper, we consider a class of DNLS equations 
featuring a general type of nonlinearity characterized by 
a nonlocal response as well as an arbitrary polynomial de- 
pendence on the intensity, that includes the case of non- 
local saturable nonlinearity. We show that, under certain 
conditions, the PNp can be made exactly zero. Quite re- 
markably, the corresponding equations are generally not 
integrable, leading to nontrivial soliton dynamics. We 
present examples of equations that conserve the total en- 
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ergy (or norm), suggesting that the predicted effects can 
be observed in specially engineered periodic structures. 

Our presentation is structured as follows. In section II, 
we present the general setup of the continuum and dis- 
crete models of the present work. In section III, we derive 
the class of PNp-free DNLS equations and extract two 
subclasses conserving different physical quantities. The 
case of cubic (Kerr) nonlinearity is considered in detail 
in section IV. Finally, in section V, we briefly summarize 
our findings and present our conclusions. 

II. SETUP 

The PNp potential is a feature of discrete models and 
does not appear for continuous and translationally invari- 
ant nonlinear Schrodingcr (NLS) equation. We demon- 
strate below that, by performing appropriate discretiza- 
tion of the continuous NLS equation, it is possible to 
obtain a broad class of DNLS models where PNp is also 
absent. We note that a conceptually similar approach 
has been developed for Klein-Gordon type lattices [20- 
27] and is now starting to emerge in NLS settings as 
well [28]. We also note in passing the interesting con- 
current efforts but in a different direction (namely that 
of obtaining exact solutions with a free parameter and 
demonstrating the absence for those of PN barriers) of 
[29 31]. Here, we present our methodology for the gen- 
eralized NLS equation of the form 

#t + ^Vxx + G'(|V|')^ = 0, (1) 

where 'ipi^t t) is a complex function of two real variables; 
G(^) is a real function of its argument and G'(^) = 
dG/d^. 

We introduce the lattice sites Xn = nh, where h is the 
lattice spacing and n — 0,±1,±2, ... We also introduce 
the following shorthand notations 

i){Xn-\) = t/j-, tp{Xn) = Ipn, tp{Xn+l) = '4'+, (2) 

and will focus only on discretizations that involve such 
nearest neighbor sites. 

Our more specific aim will be to construct the discrete 
analogues of Eq. (1) of the form: 

#„-|-r(^/'_,^*,V'„,V'n.V^-H:'0+) = 0, (3) 
such that the ansatz 

Mt) = /ne*"*, (4) 

reduces Eq. (3) to the three-point discrete problem of 
the form 

-UjU + R{f-,fn,f+) = 0, (5) 

whose solution can be found from a reduced two-point 
discrete problem u{f-,fn) = 0. Such a selection will 



entail a mono-parametric freedom for the resulting alge- 
braic problem leading to stationary state solutions. This 
will, in turn, be responsible for the effective translational 
invariance in what follows. Such models will be called 
the models free of Peierls-Nabarro potential (PNp-free 
models, for short) since their stationary solutions can be 
placed anywhere with respect to the lattice. 

Following the method proposed in [28], we derive a 
wide class of PNp-free DNLS equations. Then we demon- 
strate that they contain subclasses conserving the classi- 
cal norm 

iV = l]|V'n|^ (6) 
n 

or modified norm 

^ oo 

^=2 E (7) 

n= — oo 

and classical momcntmn 

CSO 

p = i Yl (V'nV^; - V':v'+) • (8) 

n=—oo 



III. PNP-FREE DNLS EQUATIONS 
A. Auxiliary problem 

Firstly, we formulate an auxiliary problem. Seeking 
stationary solutions of Eq. (1) in the form 

ij{x,t) = f{x)e'^\ (9) 

we reduce it to an ordinary differential equation for the 
real function f{x), 

D{x) = f" - 2cjf + 2fG'{f) = 0, (10) 

having the first integral 

u{x)={ff-2ujf + 2G{f)+C = 0, (11) 

where C is the integration constant. 

We then identify discretizations of Eq. (10) of the form 

D{f_,fn, U) = 0, (12) 

such that solutions to the three-point discrete Eq. (12) 
can be found from a reduced two-point problem 

u{f-,fn) = ^{fn- f-f 
-2a;/_/„ + 2C?(/_,/„) + C = 0, (13) 

which is a discrete version of Eq. (11), assuming that 

G{f-, fn) reduces to G{f'^) in the continuum limit {h 
0). In the present study, we will set C = 0, which is 
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sufficient for obtaining the single-humped stationary so- 
lutions. However, the case of arbitrary C enables one to 
construct all possible stationary solutions to the corre- 
sponding discrete model [32]. 

Taking into account that Eq. (10) is the static Klein- 
Gordon equation with the potential 

V{f) = wf - G{f), (14) 

a wide class of discretizations solving the auxiliary prob- 
lem has been offered in the very recent work of [23, 24] . 

For example, discretizing the left-hand side of the iden- 
tity {l/2)du/df = D{x), we obtain the discrete version 
ofEq. (10), 



£)i (/_,/„,/+) 



u{fnj+) - u{f-,fn) 
f+-f- 



0. (15) 



Formally, Di{f-, fn, f+) = is a three-point problem 
but, clearly, its solutions can be found from the two-point 
problem u{f-, /„) = and thus, the auxiliary problem is 
solved. We note, in passing, that this type of argument 
was first proposed in [22]. 

Before we come back to our main problem of finding 
PNp-frcc discretizations for Eq. (1), we should remark 
that among the solutions to the auxiliary problem we 
should select the ones which can be rewritten in terms 
of ipn and V'n in the desired form of Eq. (3). This can 
be done easily, e.g., when Di given by Eq. (15) is writ- 
ten in a non-singular form (i.e., when the denominator 
cancels with an appropriate factoring of the numerator) . 
This always occurs if G(^) is polynomial and if fn) 
possesses the symmetry 

«(/_,/„) =w(/„,/_). (16) 

We thus focus on G{^) in the form of Taylor expansion, 



G(l^l^) 



fc=i 



2k 



(17) 



with real coefficients Ofe. 

The most general polynomial, two-point discrete ver- 
sion of Eq. (17), possessing the symmetry G(/_,/„) = 
G(/„,/-)is 



G{f-Jn) 



k=l 



akQ2kif-, fn), 



(18) 



where Q2k{f-, fn) involves only the terms of order 2k: 

Q2k{j-,jn) = =^fc , (ly) 

with Qi^k being free parameters such that X]i=o 0- 
Then we substitute Eq. (13) with G{f-,fn) given by 
Eq. (18) into Eq. (15) to obtain 



-U>fn +2h2if-- 2/n + /+) 



+5E 



fe=i 



akS2k-l{f-, fn, f+) 

Ek 
i=0 1i,k 



= 0, 



(20) 



where the function 

S2k-l{f-, fn, f+) = <l0,k ^ f-f+'' ^ 
m=0 

fc 2fe-i-l 

+ E*.'^(/n E f-fl'-'-'- 
i=l m=0 

i-l 



2fe-l 



m=0 



(21) 



involves only the terms of order 2fc — 1. 



B. Main problem 

PNp-free discretization of NLS equation Eq. (1) has 
the form: 



OO 

+5E 



fe=i 



akS2k-l{'ip-,-'Pn,''P+) 

Ek 



0, 



(22) 



where S2k-i{ip-, ij^n, V-'+) is any function that, upon sub- 
stituting Eq. (4), reduces to 52fc-i(/-, /n, ./+)e*'^*, with 
S2k-i{f-, fn- f+) given by Eq. (21). Indeed, the sta- 
tionary solutions to Eq. (22), satisfying the three-point 
problem, Eqs. (20) and (21), can be found from the re- 
duced two-point problem, Eqs. (13), (18), and (19). 

To complete the construction of the PNp-free DNLS 
equation we need to obtain a suitable function 
S2k-i{4'-,''Pn,^+) reducible by the ansatz Eq. (4) 
to 52fc-i(/-,/„,/+)e*-*. Clearly, S^k-iif-, fn, f+) can 
have several counterparts <S2fe-i(V'-, V'ni V'+)- For exam- 
ple, both IV^nj^V-'- and VkV*- give (/^/_)e*'^* upon sub- 
stituting Eq. (4). 

As it can bo seen from Eq. (21), the typical term of 
S2k-i{f-,fn,f+) is ^1^^+ with k + l + m = 2k-l, 
and k,l,m> 0. This term can be transformed, e.g., to 



for A: > 0, / > 0, and m > 0, correspondingly. 
It can also be transformed, e.g., to 



(23) 



(24) 



f':fLfT^r-rn\r--'^L-'^7l 

for A; > and / > 1, or to 

for A; > 1 and I > 2, and so on. 

One can see that the number of possibilities rapidly 
increases with increase in the order of the term, 2k — 1. 
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C. PNp-free DNLS equation conserving N 



where 



Let us consider the following PNp-free DNLS equation 
1 



CXD 

It. 



fe=i 



akS2k-i{'>P-,tpn,tp+) 

Ek 



= 0, 



(26) 



where 



S2k-i{i^-, 11^71,^+) 

k 2fe-i-l 



m=0 

i-1 



+ 



V'n'"'-'|Ek-V';"'"'"|)- (27) 



m=0 



For this model, the classical norm of Eq. (6), is con- 
served. 

To construct this equation we had to set go,fc = to 
drop the first sum in the right-hand side of Eq. (21), 
because it does not contain /„. 

It is possible to construct some other DNLS equations 
conserving N and one example will be given later for the 
Kerr nonlincarity. However, the above equation is inter- 
esting because it involves only the on-site nonlinearities 
modified through inter-site coupling. 

Stationary solutions to Eq. (26) can be found from the 
two point problem Eq. (13) with G(/_, /„) given by Eq. 
(18), where 

Ek / _ri r2k—i, f2k—ifi\ 
, , , r r , — i=l'ii,k \J -Jn -rj- Jn) ,cq\ 
^2k[J-,Jn)- — , (28) 

^2^i=l <li,k 

i.e., in comparison with Eq. (19), we simply drop the 
term with i = 0. 



S2k-l{f-, fn, f+) 



k-1 



(/_+/+){go,.E/-"/f"'""'^ 

m=0 

Lfc/2J k-i-l 

E( f2i f2m f2(k-i-m-l) 

q2i,k\Jn J- J+ 

m=0 



i=l 



(31) 



m=0 



The resulting dynamical model (corresponding to (30) 
and Eq. (31)) of the form 



1 akS2k-l (V'- ,^n,^l^+) 



(32) 



with 



S2k-lilp-,'>Pn,'>P+) 

k-1 

) 



(v.-+^+){go..i:k->f"'""' 

m=0 

Lfc/2J k-i-l 
■ ' m=0 



i=l 



m=0 



(33) 



respects the conservation laws of both the modified norm 
of Eq. (7) and the momentum of Eq. (8). Stationary 
solutions to Eq. (32) can be found from the two-point 
problem Eqs. (13), (18), and (29). 



D. PNp-free DNLS equation conserving N and P 

In Eq. (19), let us drop all the terms corresponding to 
odd i . We can write the result in the form 

Lfc/2J „ ( f2i f 2fe-2i I r2k-2i f2i\ 

(29) 

where [A;/2J is the largest integer not greater than fc/2. 
Substituting Eqs. (13), (18), and (29) into Eq. (15) we 
obtain 



-^/n +2}^(f-- 2/n + /+) 
, 1 akS2k-l{f-, fn, f+) _ p, 

+ 2 ^ ^lk/2i - ^"^"^ 

fe=l Z^i=0 y2s,fc 



IV. CUBIC (KERR) NONLINEARITY 

A. Examples of PNp-free models 

Let us study in detail Eq. (1) with Kerr nonlinearity, 
i.e., with G'(IVP) - IV'P, 



(34) 



We then have G(|V'P) = |V'lV2, i-e., in Eq. (17), a2 = 
1/2 and afc = for fc ^ 2. 

To construct PNp-free DNLS equations we write Eq. 
(20) and Eq. (21) for the case of Kerr nonlinearity as 

-./„ + 5ij(/_-2/„ + A) + |^^^=0, (35) 
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and 

Ssif-, fn, /+) = a{fl + f-fl + f_f+ + f_) 

+27(/+ + /-)/^, (36) 

where we introduced the foUowing shorter notations for 
the free parameters: a = go, 2, P = (?i,2! and 7 = q'2,2- 

Solutions to the three-point problem of Eq. (35) can be 
found from the following two-point problem [Eqs. (13), 
(18), and (19)] 

1 (/„ - f_f - 2a;/_/„ 

, a(/l + ft) + P{f-fl + f-fn) + 2^f_fl _ ^ 
+ 2(a + /? + 7) " ^ ^ 

The ensuing PNp-free DNLS equation with Kerr non- 
linearity is 



where iSa (■)/'-, V'n; "0+) is any function that, upon sub- 
stituting Eq. (4), reduces to e^'^'S'a (/_,/„, /+) [i.e., 
respecting the phase invariancc of the equation], with 
Sz{J-,fn,f+) given by Eq. (36). Some guidelines on 
how to construct Ss{ijj-,tjjn,tlJ+) can be found in Sec. 
IIIB. 

From Eq. (26) and Eq. (27), at 03 = 1/2 and Uk = 
for fc ^ 2, we obtain the Kerr-type DNLS equation 
conserving the classical norm N 



+ 



4(/J + 7) 



ii^n + ^ ~ ^V'n + V'+) 

■/3(|V-P + |^-V+l + IV'+P + IV'nP) 

+27(|V-Vn| + |^^nV'+l)l =0. (39) 



The two-point equation for finding the amplitudes of sta- 
tionary solutions to Eq. (39) is Eq. (37) with a = 0. 

Similarly, from Eq. (32) and Eq. (33), we obtain the 
Kerr-type DNLS equation conserving modified norm N 
and momentum P 



ii^n + ^ (V- 



+ 



4(q! + 7) 



a(|^_P + |V;+|2)+27|VnP =0. (40) 



Amplitudes of stationary solutions of Eq. (40) satisfy Eq. 
(37) with j3 = 0. 

Notice that the intcgrable discretization of [15, 16] is 
obtained from Eq. (40) as the special case of a = 0. 
For a ^ 0, this model can be regarded as a Salerno- 
type model [33], i.e., a homotopic continuation including 
the intcgrable limit and reducing to NLS equation in the 
continuum limit. 



DNLS equations of the form of (39) and (40) do not, of 
course, exhaust the list of possible PNp-free models with 
Kerr nonlinearity. To give one more example, we note 
that the last term of Eq. (36) can be used to produce 
the following DNLS equation conserving classical norm 

ii'n + ^ (V'- - 2V'n + ^+) + \ (V'- + i>+) IV-nj' 

+l{r-+r+)i>i = o- (41) 

Amplitudes of stationary solutions to Eq. (41) can be 
found from Eq. (37) at a = (3 = 0. 



B. Soliton solutions 

We now compare some properties of the classical DNLS 
equation, 

model I : 

itpn + ^(V'- - 2V'n + V'+) + IV'nPV'n = 0, (42) 

with these of the Af-conserving model of Eq. (39) with 
(3 = 0, 

model // : 



+ ^(|V-Vn| + l^^nV'+l)=0; 



(43) 



and those of the N- and P-conserving model of Eq. (40) 
with 7 = 0, 

model /// : 



+ 



ii'n + ^ ~ 2V'n + 

^^(i^_p+i^,n=o. 



(44) 



All three models share the same continuum limit, the 
intcgrable NLS equation Eq. (34), and thus, in the 
regime of weak discreteness (small lattice spacing h), 
their soliton solutions of the form of Eq. (4) can be ex- 
pressed approximately as 



cosh.[Ah{n — xo)] 



eM-i{AV2)t], (45) 



where A and u> = are the soliton amplitude and 

frequency, respectively. 

The approximate solution of Eq. (45) contains the free 
parameter Xq defining the soliton position. However, in 
contrast to the NLS equation of Eq. (34), where xq can 
be chosen arbitrarily due to translational invariance, the 
DNLS models usually have stationary soliton solutions 
only for a discrete set of values of a;o (e.g. on-site, xo = 0, 
and inter-site, xq = 1/2). This is true, for example, for 
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the classical DNLS of model I and for the Salerno model 
[33], among others. The models // and ///, by construc- 
tion, are among the members of a wider class of DNLS 
equations proposed in this paper, where stationary soli- 
ton solutions exist for any xq, or, in other words, they 
can be placed anywhere with respect to the lattice; oth- 
erwise put, the Peierls-Nabarro potential is absent for 
stationary solutions of these models. 

Let us now describe the exact soliton solutions to the 
models /, //, and ///. 

An explicit formula does not exist for the stationary 
soliton solutions of model /. Such solutions can be ob- 
tained using the fixed point algorithms [34], but, as men- 
tioned above, only for xq = and xq = 1/2. 

Model II, model of Eq. (41), and also model of Eq. 
(40) at a = 0, have the same equations for the ampli- 
tudes of stationary solutions. However, the latter model, 
as mentioned above, is the integrable DNLS equation 
[15, 16] and thus, an exact stationary solutions for these 
models can be obtained explicitly in the form 



1 



sinh 



h cosh[/i(n — xo)] 



exp 



(46) 



where xq is the parameter defining the soliton position 
and it can obtain any value from [0,1). The soliton 
frequency w = /i~^(l — cosh/z) and amplitude A = 
h~'^ sinh^ H are expressed in terms of the free parameter 
M > 0. 

Model /// has the solutions of the form of Eq. (4) 
with /„ derivable from the two-point problem 



1 



(/„-/_)' -2a;/_/„ + ^ (/I +/^) = 0. (47) 



The soliton can be constructed by setting an arbitrary 
value for /_ (or /„) in the range [Am, As] and finding /„ 
(or /_) from the quartic Eq. (47). Quantities A^ and 
As are the amplitudes of solitons centered between two 
lattice sites and on a lattice site, respectively. We have 
Am — and As can be foimd from the condition 

that two distinct real roots of Eq. (47) merge into a 
multiple root. The arbitrariness in the choice of initial 
value of /_ (or /„) implies the absence of the Peierls- 
Nabarro potential and the possibility to place the soliton 
anywhere with respect to the lattice. 



C. Soliton's internal modes 

Let us stiidy the stability of stationary soliton solutions 
for the models described in Sec. IV B. In this study we 
calculate the soliton's internal modes and frequencies of 
these modes. 

Following the methodology of the paper [35], to study 
the stability of the solution Eq. (4), we consider the 
complex perturbation e„(t) in the frame rotating with 
the periodic solution: 



Substituting Eq. (48) into the classical DNLS equation 
Eq. (42) (model /) we find that the linearized equation 
satisfied by e„(f) is 

tin -^^n + ^ (e- - 2e„ -h e+) + 2f^en + f^e^ = 0. (49) 

Similarly we obtain the linearized equations for e„(t) 
for the A/'-conserving model II of Eq. (43), 



1 

2^2' 



ten - we„ -I- (e_ - 2e„ -|- e+) 



+ ^fn{f-+f+)[en + M^n)] 

+ i/2[Re(e_)+Re(e+)] = 0, (50) 

and for the N- and P-conserving model /// of Eq. (44), 

iin - uJCn + ^(e- - 2e„ -I- e+) 
+ i(/_ + /+) [/-Re(e_) + /+Re(e+)] 

+ \{f_ + fl){e.+e+)=0. (51) 



tions can be written as 



where 



Mt) = [fn + en{t)] 



(48) 



Defining e„(f) = a„(t) + ibn{t), the linearized equa- 

a\_/o n \ ( & 

vectors a and b contain a„ and 6„, respectively. Sta- 
tionary soliton solutions are linearly stable if and only if 
the eigenvalue problem, det(Jfl — A^I), has only real and 
nonpositivc solutions for [35]. 

In our stability analysis, a soliton, having a frequency 
w = 1, is placed at the middle of chain of 400 sites with 
periodic boundary conditions. Different magnitudes of 
the discreteness parameter h and different positions of 
solitons with respect to the lattice, xq, are considered. 

Most of eigenfrequencies V— A^ appear within the band 
between u) and w + where w is the frequency of 

soliton and is the maximum frequency of the lin- 

ear spectrum of trivial solution, /„ = 0. We do not 
show these eigenvalues in the figures. Eigenvalues ap- 
pearing outside of this band are related to soliton's in- 
ternal modes. The spectra of solitons in the models /, 
//, and /// always contain two zeroes. However, spectra 
of solitons in the PNp-free models (including models II 
and ///), due to their effective translational invariance, 
always contain two additional zeroes. 

Results for the classical DNLS equation (model /) are 
presented in Fig. 1 for (a) xq — Q and (b) xq — \/2. In 
(a), the on-site soliton is stable because all A^ are real and 
nonpositive, while in (b), the inter-site soliton is unstable 
because we have positive A^ (shown by open circles). The 
solid line shows the bottom edge of the linear band and 
dotted line shows the always existing eigenvalues A^ = 0. 

In Fig. 2, the results for the A-conserving, PNp-free 
DNLS equation (model //) are presented. The always 
existing eigenvalues A^ = are shown by dotted line and 



7 



7 0.6 



(a) 




0.4 0.8 1.2 1.6 2.0 



0.4 0.8 1.2 1.6 2.0 




FIG. 1: Model I, spectrum of the soliton with frequency uj = 1 
(bottom edge of the continuous spectrum) for (a) stable on- 
site and (b) unstable inter-sitc configurations. Spectrum al- 
ways contains a pair of zero- frequency modes shown by dashed 
lines. For h < 0.4 discreteness is weak and the soliton in 
(a) has a mode with nearly zero frequency corresponding to 
the translational mode of the continuum NLS equation. The 
frequency of this mode grows rapidly with increase in h for 
h > 0.4 and it enters the continuous band at /i = 0.75 where 
soliton loses its mobility. 



the additional two zeroes are shown by dots; the latter 
reflect the translational invariance of the soliton. Note 
that panel (b) corresponds to the soliton having a;o = 
1/4, i.e., placed non-symmetrically with respect to the 
lattice points. For h > 1.25, below the linear spectrum, 
there exists a soliton internal mode. This mode was not 
observed for the soliton placed at xq = 1/2. 

Similar results for the N- and P-conserving, PNp-free 
DNLS equation (model ///) are shown in Fig. 3. In 
contrast to models / and II, soliton in model /// has 
the internal modes lying not only below, but also above 
the linear band. These modes are shown in Fig. 4 by 
dots and the upper edge of the linear spectrum is shown 
by the solid line. 



FIG. 2: Model //, spectrum of the soliton with frequency 
uj — 1 for (a) on-site and (b) asymmetric configurations. Both 
configurations are stable. A soliton internal mode bifurcates 
from the bottom edge of the continuous frequency band at 
rather large h. In addition to an always existing pair of zero- 
frequency modes (shown by dashed lines), for any xo, there 
exist another pair of zero- frequency modes (shown by dots), 
rcficcting the translational invariance of the soliton and the 
absence of the Peierls-Nabarro potential. 




D. Mobility of solitons 

Solving the eigenvalue problem formulated in Sec. 
IV C, we find solutions of corresponding DNLS equa- 
tion of the form of Eq. (48). Solutions are accurate 
for small amplitudes of the eigenvectors, e„(t). Here we 
examine the solutions corresponding to the translational 
eigenmode (since "kicks" along this eigendirection may 
be responsible for/related to motion in these lattices). 

We define the position of discrete soliton as its center 
of mass: 



FIG. 3: Model ///, bottom part of the spectrum of the soli- 
ton with frequency cu = 1 for (a) on-site and (b) inter-site 
configurations. Both configurations arc stable. Soliton inter- 
nal mode bifurcates from the bottom edge of the continuous 
frequency band. For any xo, the spectrum contains two pairs 
of zero-frequency modes, which is a common feature for the 
PNp-free models. 

where Nn = \tpn\^', then the soliton velocity is u = dS/dt. 
This definition is used for TV-conserving models / and 
//, and for A^-conserving model ///, in Eq. (52), we 

naturally use = (l/4)(V'„V'-+V'^'0-+V'nV'+ + V'n'^-i-) 
instead of Nn- 

Setting initial conditions according to Eq. (48) with 
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FIG. 4: Same as in Fig. 3 but for the top part of the spectrum. 
In contrast to models I and //, the soliton in the model /// 
has internal modes with frequencies lying above the continu- 
ous band. Solid line shows the upper edge of the continuous 
band. 
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FIG. 5: Model I, h = 0.4. Time evolution of (a) soliton 
velocity, (b) norm, (c) modified norm, and (d) momentum. 
Soliton frequency is w = 1. 



sufBciently small amplitude of translational eigenmode 
e„(t), we then integrate numerically the DNLS equations 
to study the soliton mobility at different h. 

For all three models we boost the solitons initially 
placed at xo = 0. 

As Fig. 1(a) suggests, the soliton of model / has a 
translational mode with nearly zero frequency only for 
h < 0.4. For such a small h, discreteness is weak and 
model / can be regarded as weakly perturbed continuum 
NLS equation, which supports moving solitons. However, 
the frequency of the translational mode increases rapidly 
for ft, > 0.4 and it enters the continuum frequency band 
at ft = 0.75 when the soliton completely loses its mobility. 
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FIG. 6: Model /, h — 0.4. Time evolution of soliton velocity 
for different initial velocities. Soliton frequency is w = 1. 
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FIG. 7: Model //, h = 2.0. Time evolution of (a) soliton 
velocity, (b) norm, (c) modified norm, and (d) momentum. 
Soliton frequency is a; = 1. 
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FIG. 8: Same as in Fig. 7 but for model 



9 



In Fig. 5 we show the time evolution of sohton's veloc- 
ity, V, norm, N , modified norm, N, and momentum, P, 
for the soliton of frequency w = 1 in model I at h = 0.4. 
Soliton was boosted with the initial velocity ofv = 0.024. 
One can see that norm is conserved while all the other 
parameters oscillate while the soliton propagates along 
the chain. The soliton's velocity is maximum (minimum) 
when it passes the on-site (inter-site) configuration with 
minimum (maximum) energy. Its average speed is equal 
to 0.016. 

To boost the soliton in model I one has to apply suffi- 
ciently large initial momentum (i.e., larger than a thresh- 
old in order) to overcome the Peierls-Nabarro potential. 
This feature is illustrated by Fig. 6, where we boost the 
soliton with different initial velocities: (a) vq = 0.012, 
(b) vo = 0.018, and (c) vq = 0.024 (same as in Fig. 5). 
The soliton does not propagate in (a) and in (b), instead, 
it oscillates near the minimum energy configuration. 

The soliton kinematics/dynamics in the PNp-free mod- 
els // and /// are qualitatively different from those in 
model / because soliton is not trapped by the lattice and 
thus can be accelerated by even weak external fields. 

In Fig. 7 and Fig. 8 we show the results for models II 
and ///, respectively, similar to that presented in Fig. 
5 for model /. Note that for model / we used a rather 
small discreteness parameter ft, = 0.4 while solitons in 
Fig. 7 and Fig. 8 propagate along highly-discrete chains 
at /i = 2. It can be clearly seen that model // conserves 
classical norm N, while model /// conserves modified 
norm N and momentum P. 



V. CONCLUSIONS 

We have described a general and systematic method of 
constructing spatial discretizations of NLS-type models, 
whose stationary soliton solutions can be obtained from a 
two-point difference problem. In this setting, finding sta- 
tionary solutions becomes tantamount to solving simple 



nonlinear algebraic equations. We have also illustrated 
the connections of the resulting models with the inte- 
grable discretization of the NLS equation, of which they 
are a natural generalization for cubic nonlinearities (our 
construction was given for arbitrary polynomial nonlin- 
earities of a particular parity); furthermore, the differ- 
ences of such models from the standard discretization of 
the NLS equation often encountered in physical applica- 
tions have been highlighted, both in terms of the relevant 
dynamical (solitonic) behavior as well as in terms of the 
underlying conservation laws present in the various mod- 
els. 

For the case of cubic nonlinearity we demonstrate that 
in the constructed PNp-free DNLS chains solitons are 
stable in a wide range of the discreteness parameter h. 
Semi-analytical moving soliton solutions are found in the 
regime of propagation at slow speed. Several numerical 
examples illustrate the qualitative difference in slow soli- 
ton dynamics in the classical DNLS equation and the 
PNp-free DNLS equations. In the classical model to 
boost the soliton one has to apply sufficiently large initial 
momentum to overcome the Peierls-Nabarro potential. 
In the PNp-free model, the soliton behaves differently 
and can be accelerated by even weak external fields. 

We believe that our results may suggest novel possi- 
bilities for engineering nonlinear lattices optimized for 
efficient control over the localization and propagation of 
discrete solitons in various physical contexts. 
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